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Lattice studies of the infrared regime of gauge theories are complicated by the required extensive 
limits, the performed gauge fixing and the demand for high statistics. Using a general power 
counting scheme for the infrared limit of Landau gauge SU(N) Yang-Mills theory in arbitrary di- 
mensions we show that the infrared behavior of Greens functions is both qualitatively and quanti- 
tatively similar in two, three and four spacetime dimensions. Therefore, lower dimensional lattice 
simulations can serve as a simplified laboratory to analyze the presently applied approximations 
and to obtain first results for higher correlation functions. 
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Introduction: Lattice studies of the infrared (IR) regime of Yang-Mills theory, see e.g. [jlj], are 
important to understand its genuinely non-pertubative aspects. This is particularly promising since 
the detailed mechanism for gluon confinement could be connected to directly measurable correla- 
tion functions of the colored fields within the Kugo-Ojima [||] and Gribov-Zwanziger ^ scenarios. 
These rely on a strong IR increase of the ghost dressing function that was predicted by continuum 
methods 0, [5[ 0] and has been clearly confirmed in lattice simulations ^ . Combined efforts via 
all available methods such as Dyson-Schwinger equations (DSE) [Q, ||, ^ |7|], renormalization group 
(RG) techniques [[|] and lattice gauge theory studies [jl]] led during the last years to a coherent pic- 
ture of the IR regime in Landau gauge. The continuum approaches allow in particular to analyze 
the IR scaling limit of general Greens functions [[J. These are compatible with a dominance of 
the gauge fixing part of the action, as argued for in JTQp . A simultaneous analysis of DSE and RG 



methods allowed furthermore to show that this IR fixpoint is unique JllQ. 

Recently, challenging lattice results by several groups on larger lattice sizes [12] suggested that the 
IR scaling of the gluon propagator might be weaker than predicted by continuum methods and does 
not show the IR vanishing required within the Gribov-Zwanziger scenario. The finite size correc- 



tions have been previously estimated within a Dyson-Schwinger study on a compact manifold [ |13| ] 
and suggest that despite significant corrections the expected IR behavior should be visible with the 
quite remarkable lattice sizes employed. However, the four dimensional studies still feature signif- 



icant statistical fluctuations in the IR which are less pronounced in two [14] and three dimensions 



[15] where deviations from the continuum predictions are not observed. 

The discrepancy of the new lattice data with the continuum results can basically have two reasons. 
First, it is possible that the IR regime requires even more refined numerical methods on the lattice 
side. As we will argue below, in this case lower dimensional studies should allow to test the in- 
volved approximations in a substantially simpler setting. This conclusion stems from a study of the 
dependence of the IR limit on the spacetime dimension d. We establish a manifest power counting 
scheme for general vertex functions of Landau gauge Yang-Mills theory in arbitrary dimensions 
and find qualitatively similar results as in the four-dimensional analysis given in [g]. 
The other reason for the discrepancy could be that the IR behavior predicted by continuum meth- 
ods - and in particular the value of the IR scaling parameter fc - suffers from the presently studied 
approximations. In an approximation based on the propagator DSEs and using a bare ghost-gluon 
vertex it takes a value K w 0.595 in four dimensions. A possible dressing of this vertex could 
change this value [Q], and the continuum prediction would be compatible with the recent lattice 
results if the IR exponent k is slightly smaller. Therefore, we also discuss the dependence of the 
dynamical building blocks in the Dyson-Schwinger equations on K. 

IR exponents for arbitrary d: In the following we will perform a scaling analysis for the IR 
regime of Yang-Mills theory in arbitrary dimensions. In contrast to the generic IR limit we study 
here the limit where the coupling g is kept fixed and has no inherent scaling dependence, which is 



more directly accessible in lattice simulations and has been analyzed in Q14J , [15J. As in four dimen- 
sions [9], the starting point for the IR analysis is the non-renormalization of the ghost-gluon vertex 
in Landau gauge [|I6|], which implies a finite vertex in the IR [0J. This property depends purely on 
the transversality of the gluon propagator and is therefore valid in arbitrary dimensions. The IR 
behavior of the propagators and vertices well below its inherent scale g 2 l ( 4 ~ rf ) (respectively Aqcd 
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Figure 1: The DSE for the ghost propagator. 



in d = 4) is determined via renormalization group arguments by scaling relations. The propagators 
of the gluons and ghosts 
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are given in terms of dressing functions whose IR behavior is described by a power law ansatz 

Z(p 2 )=c Q ,2-(p 2 ) So2 , G{p 2 )=c 2fl -{p 2 t-\ (2) 

and similar for the vertices. Here we denote the IR exponent of a vertex with In ghost and m gluon 
legs by &in, m and the corresponding coefficient by c 2 „, m . Whereas this coefficient is a constant for 
the propagators it is generally a function of 2n+m— 1 momentum ratios. 

For the integral on the right hand side of the ghost propagator Dyson-Schwinger equation, cf. fig. 
H| one can use the standard expression [|6[ |n| ] 
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which shows that it scales proportional to the external momentum. The left hand side of the ghost 
DSE, which consists only of the inverse dressed ghost propagator, is proportional to (/? 2 )~^ ,0+1 , 
where the 1 is the canonical dimension. The scaling dimensions on the right hand side are d/2 from 
the integral, 5o,2 — 1 from the gluon propagator, 82 ,0 — 1 from the ghost propagator, 1/2 from the 
bare ghost-gluon vertex, and 1/2 from the dressed ghost-gluon vertex which features no anomalous 
scaling. Defining the parameter K as K := — 5 2j o> this condition on the IR exponents yields in four 
dimensions the well-known result 5q, 2 = 2k. Ind dimensions we have [|5[ ® 



§2.1 



30,2 



2K + 2--. 



(4) 



In order to transform the infinite hierarchy of Dyson-Schwinger equations into a closed system we 
perform a skeleton expansion. This yields an infinite tower of graphs involving only primitively 
divergent vertices which we will analyze in a first step. The first order of the skeleton expansion of 
the DSE for the three-gluon vertex is depicted in fig. |2[ 




Figure 2: 1-loop part of the skeleton expansion of the DSE for the three-gluon vertex. 
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We start with the ghost triangle which turns out to be one of the IR leading diagrams. Subtracting 
the canonical dimension 1/2 we get from the IR counting rules the anomalous IR exponent: 



8^= d - + ?>{-K-l) + 3\- l - = -3K+ d --2. (5) 

For the four-gluon vertex we can apply the same procedure. In fig. |3| we show the first order of its 
skeleton expansion. Here, again the ghost rectangle gives one of the IR dominant contributions: 

8lf = d - + A{-K-\)+A\ = -AK+ d --2. (6) 

With these results for the scaling of the primitively divergent vertices it is easy to check that the 
other diagrams in the DSEs fig. || and [3] have either the same scaling or are subleading in the 
IR. Similarly one can analyze the IR exponent of an arbitrary rc-point function. This is done by 
counting the IR exponents of the individual building blocks that arise in a general loop correction. 
The resulting expressions can be simplified via topological identities. It is then straightforward to 



show that the leading IR exponent of an « -point function is given by the simple expression [18] 



d 2 n,m = (n-m)K+(l-n)( d ;-2\ . (7) 

This relation verifies a posteriori the assumption discussed above fig. [I]: the ghost-gluon scattering 
kernel is exactly as IR divergent as to make the ghost-gluon vertex finite in the IR, and this is 
independent of the value of the dimension d. 

It remains to discuss higher order corrections in the skeleton expansion which can be generated via 
insertions into given diagrams. The change of the IR exponent due to an exemplary insertion is 
given by: 

X x 

^ AAAA 4 (2k+1_!) + 2(-k-1) + 2£ + !=0, 
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i.e. the IR-scaling is unchanged. This can also be shown for all other possible extensions [|TJ]. Cor- 
respondingly, the skeleton expansion works independent of the number of spacetime dimensions 
d [Q, HI ||, [| The same holds for the IR limit of the momentum dependent couplings defined 
through the different vertices: they only feature a scaling via their canonical dimension. 
Finally, let us compare the IR behavior of the Yang-Mills Greens functions in different dimensions. 
We show the IR exponents of the ghost and gluon propagators and the three- and four-gluon ver- 
tices in table []} 



Dimension 




4 


3 


2 


K 




0.5953... « 0.6 


0.3976... 0.4 


0.2 


Ghost 


-K--1 


-1.6 


-1.4 


-1.2 


Gluon 


2K+l-d/2 


0.2 


0.3 


0.4 


3-gluon 


-3K + d/2-3/2 


-1.3 


-1.2 


-1.1 


4-gluon 


-Ax + d/2-2 


-2.4 


-2.1 


-1.8 



Table 1: The dimension dependence of the IR behavior of Yang-Mills Green functions. For each d only the 
continuum solution for K is displayed that has been qualitatively confirmed by lattice simulations. 



Whereas the situation in four dimensions is not fully conclusive, yet, the other solutions have been 
found in lattice studies in two and three dimensions. These results for the IR exponents of the prop- 



agators, the ghost-gluon vertex [ 19] and the three-gluon vertex in two dimensions [ 14] agree within 
errors with the corresponding values given in table [l[ The qualitative behavior does not change in 
different dimensions and even the quantitative values are very similar. Apart from the listed values 
of K there is a second branch of solutions of the DSE equation that has not been seen in lattice 
simulations. We will give some additional arguments below that this branch may be unphysical. 

Dependence of the loop integrals on the IR exponent: The DSEs involve loop integrals over 
bare and dressed vertex functions. As discussed the dressing functions exhibit a power law with 
appropriate IR exponents. The corresponding coefficients depend on the IR exponent K. For the 
respective coefficients of the propagators, eq. (|2|), which can be computed from the 2-point in- 
tegral, eq. (||), this (("-dependence is shown in fig. |J where ^-independent prefactors have been 
dropped. The interesting point is, that although the curves in fig. || differ considerable, the physical 
values obtained from the Dyson-Schwinger solution are far away from the poles and thereby the 
^-dependence in their vicinity is rather mild and qualitatively similar in each case. 
In general the vertex integrals feature more complicated tensor structures with an increasing num- 
ber of independent tensor components. However, there are general methods to decompose such 
tensor integrals to standard scalar integrals [|20|]. The tensor integrals in the DSE for the ghost- 
gluon and three-gluon vertex reduce to 3 -point integrals of the form 

r d d k 1 1 1 

h(P,q)=J ^ 1%)d ^ k + p yy l ((k-q) 2 ) v 2 (£2)v 3 ^ 

A general expression for such scalar one-loop 3-point integrals in arbitrary dimensions and with 



arbitrary powers of the propagators has been obtained in [17, 21]. An analytic expression can be 



obtained via appropriate series representations for the resulting generalized hypergeometric func- 



tions [22, 23]. Here we will discuss the dependence of the overlap of the tree-level tensor with 



5 



Lower dimensional Yang-Mills theory as a laboratory to study the infrared regime 



Kai Schwenzer 




K K 

Figure 4: The K"-dependence of the unrenormalized dressing integrals appearing in the DSEs for the ghost 
(left) and gluon propagator (right). The dotted, dashed and solid lines show the curves for d — 2,3 and 4 
respectively, see also ref. [pfl. The full and open dots represent the solutions of the first and second branch. 



the IR-dominant ghost triangle which presents the IR-leading contribution to the Dyson Schwinger 
equation for the three-gluon vertex in fig. ||. Since this loop involves only ghost-gluon vertices 
that remain bare to leading order, its IR behavior can be analyzed semi-perturbatively using the 
scaling form of the dressed ghost propagators, see eq. @. We consider the special kinematic 
configuration given by the symmetric point p 2 = q 2 = r 2 . At this point the dressing function 
of the leading ghost loop correction to the three-gluon integral scales according to eq. (|5|) as 
Zo,3 = co^(k) ■ (p 2 )- 3,c +d/2-2 w h ere ^e dependence of the coefficient on the scaling parameter K 
has been made explicit. The coefficient is shown as a function of K in fig. |5J, where the full and 
open points represent the values for the first and second branch and fc-independent factors have 
been dropped again. 



Figure 5: The fc-dependence of the 
overlap of the unrenormalized ghost 
triangle with the tree-level tensor 
at the symmetric point. The dot- 
ted, dashed and solid lines show the 
curves for d = 2,3 and 4 respec- 
tively. The full and open dots repre- 
sent the solutions of the first and sec- 
ond branch. 

The vertex integral has zeros at fc = in two, fc = 0.5 in three and K = 1 in four dimensions which 
coincide precisely with the second branch of solutions of the DSEs that is not seen in lattice simu- 
lations. Such zeros in the vertex dressing function would lead to an IR-vanishing coupling. Since 
Yang-Mills theory is apparently a strongly interacting theory, this indicates that these solutions 
might not be physically relevant. As in the case of the propagators the fc-dependence in the vicinity 
of the solutions of the first branch is similar in all dimensions, ynamics would dominate even more. 

Conclusions: As in the 4-dimensional case, the IR behavior of Greens functions can be extracted 
via a skeleton expansion in arbitrary dimensions. The IR limit of Greens functions is surprisingly 
insensitive on the spacetime dimension. As a consequence Yang-Mills theory in lower dimensions 




6 



Lower dimensional Yang-Mills theory as a laboratory to study the infrared regime 



Kai Schwenzer 



should be qualitatively similar to the 4-dimensional theory, and the confinement mechanism might 
possess identical features. Therefore, corresponding lattice simulations should provide interesting 
qualitative information for the physical case. Such studies in lower dimensions could in particular 
provide important lattice results on the IR behavior of vertex functions. This is particularly chal- 
lenging due to the observation of additional kinematic singularities in current DSE studies []22|]. 
The results on the mild fc-dependence of the remaining DSE solutions suggest that ghost domi- 
nance is a rather robust mechanism and should not depend on the details of the employed trun- 
cation scheme. This is further substantiated by the fact that the main non-linearities in the DSE 
system, which enable the non-trivial nxpoint, arise in the propagator equations whereas there is no 
non-linear feedback to infrared leading order in the vertex equations. 
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